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Abstract 



^ I Asymptotic expansions are given for large values of n of the general- 

ized Bessel polynomials The analysis is based on integrals that 

follow from the generating functions of the polynomials. A new simple 
■ expansion is given that is valid outside a compact neighborhood of the 

origin in the z— plane. New forms of expansions in terms of elementary 
functions valid in sectors not containing the turning points z — ±i/n 
are derived, and a new expansion in terms of modified Bessel functions 
is given. Earlier asymptotic expansions of the generalized Bessel poly- 
nomials by Wong and Zhang (1997) and Dunster (2001) are discussed. 
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1 Introduction 

Generalized Bessel polynomials of degree n, complex order fi and complex 
argument z, denoted by Yn{z), have been introduced in and can be 
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defined by their generating function. We have [3]: 



\ V /V n=0 

with special values 1^(0) = 1, n = 0, 1, 2, . . .. 

The generalized Bessel polynomials are important in certain problems 
of mathematical physics. For a historical survey and discussion of many 
interesting properties, we refer to [3]. 

In [2] and [13] detailed contributions on asymptotic approximations are 
given on the generalized Bessel polynomials. Wong and Zhang use integral 
representations and Dunster's approach is based on a difi^erential equation. 
Our approach also uses integrals as starting point. 

Our approach is different from that of Wong and Zhang ^13]. We give 
expansions that are similar to those of the modified Bessel functions, and the 
expansions reduce to these expansions when = 0. For the expansions in 
terms of elementary functions we give a simple description for the domains 
of validity. 

In ^we start with a simple expansion valid outside a fixed neighborhood 
of the origin. In f|3]we give expansions that hold uniformly with respect to 
z inside the sectors |ph ± z\ < ^tt — 6, where is a small positive number. 
At the end of this section we compare our results with those in |2] and [13]. 
For complementary sectors (and extensions of these) we give an expansion 
in terms of the modified Bessel function K , i(z) and its derivative. In fact, 

this is an expansion in terms of Y^{z), and this expansion holds for all z. 
In the Appendix ^ we summarize the expansions of the modified Bessel 
functions. These expansions play a role when comparing the results for the 
generalized Bessel polynomials with those for /x = 0. 



2 Representations and relations with Bessel and 
Kummer functions 



For /i = the generalized Bessel polynomials become well-known polynomi- 
als that occur in representations of Bessel functions of fractional order. We 
have in terms of the modified iiT— Bessel function [Qrj: 




http://dlmf.nist.gov/10.49.E12 
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where (p)k is the Pochhammer symbol defined by 

(p)o = l, (p)fc = ^^^±^, k>l. (2.2) 

For Yn{z) an exphcit formula reads [5] 

" / \ k 
^-(^) = E U + + (l^) • (2-3) 

The simple integral representation 

(^) = T^,^ , .. , / i""^^ (l + h*) e-'dt (2.4) 



r(n + ^ + 1) Jo 

gives the representation in (|2.3|) by expanding (1 + ^zt)" in powers of z. 

Prom the generating function in (jl.ip we have the Cauchy integral rep- 
resentation 



where C is a circle with radius smaller than l/\2z\. 

When z = all polynomials reduce to unity, the first few polynomials 
being 

Yo^{z) = l, y^'^(z) = l + l(/i + 2)z, 



Y^^{z) = 1 + (^ + 3)z + i(^ + 3)(/i + 4)^2. 



(2.6) 



More values can be obtained from the recurrence relation 

AnY^^,{z) = BnY^^,{z) + CnY>^{z), (2.7) 

where 

An = 2(2n + fi + 2){n + ^l + 2), 

Bn = {2n + n + 3)(2/i + z{2n + /i + 4)(2n + /u + 2)), (2.8) 
Cn = 2(n + l)(2n + /u + 4). 
There is also a recursion with respect to fi: 

(n + + 2)Y>^+\z) = (2n + ^ + 2 - 2/z)YJ^+H^) + (2A)lT(^), (2.9) 



and for the derivative we have 

= \n{n + ^,+ l)Y::^^{z) = !1±A±1 (yM+i(^) _ . (2.10) 

These relations are special cases of known analytic continuation formula for 
the Kummer functions. See, for example, [7J and [1(J|, p. 13]. They follow 
also from (|2.4p by integrating by parts. 

The relations with the Kummer functions are 



Y^{l/z) = {2zY+>'+^U{n + + 1, 2n + ^ + 2, 2z), 
r(n + ^ + 1) 



(2.11) 



For 5iz < it is convenient to have the representation 

y,^(-l/z) = Fi^{l/z) + U^{l/z\ (2.12) 

where 



^" <'/^* = r(2u + , + 2) + 1^ 2" + " + 2^ 

^n^(lA) = 7 + 1, 2n + ^ + 2, 2z). 

r(n + /i + 1) 



(2.13) 



For /Li = we have 

F0(z) = ^/2^e-%+i(l/z), t/0(^) = (-l)",/^e-^K„^.i(z), (2.14) 

2 V VT 2 

and this corresponds to the relation 

Yni-z) = J^e-'/^ U-irK^^M/z)+7rI^^^{l/z)) . (2.15) 

3 A simple expansion 

The interesting region in the z plane for uniform asymptotic expansions is 
a neighborhood of size 0{l/n) of the origin, where the zeros appear. For z 
outside a fixed neighborhood a simple expansion will be derived. 
First we mention 



yM(^) = 2^+l (^^y e^/^ 



X 



24n 
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This is derived in [Ij and mentioned in [3l p. 124] and [13]. More terms in 
this expansion can be obtained, for example by using the Cauchy integral 
given in (j2.5p . 

In this section we derive a simple asymptotic expansion related to the 
result in (j3.ip by expanding part of the integrand in (j2.5p in powers of 
W = \/l — 2zw. First we notice that the main asymptotic contributions 
from the contour integral in (12. 5p come from the singular point w = 1/(22:), 
and when w ~ 1/(2^) the quantity W is small. 

We have 

( V ^2n,/ii+Vi=2l^) = 2^ei/^(l + Wr^'e-^'' (3.2) 

and we expand for m = 0, 1, 2, . . . 

m—1 

(1 + Wy^'e-^/' = ^ Ll^-''{l/z)W^ + W'^U^iW). (3.3) 

fc=0 

The appearance of the Laguerre polynomials becomes clear when expanding 
both the exponential and binomial and by comparing the coefficients with 
the representation 



m=0 

Introducing this expansion in (|2.5p we find 

m—1 



Y^{z) = n! If^e'/' ^ L.'^-'il/ z)<^'^^^ + n! i?„(n), (3.5) 



k=0 



where 



and 



$w = J- /(I - 2zw)^'-'y'^ = (izr^LJhn^ (3.6) 



i?m(n) = ^ / (1 - 2u;z)(-i)/2[/„(IF) (3.7) 



with C a circle with radius less than l/|2z|. 

After the change of variable w = {\ — t/n)/{2z) we have 



f ^Ufo , t^""^^/^Um(Vt/^) 7 ^T-^- (3.8) 
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The function Um{W) is analytic in \W\ < 1 and 0(1) asW^O. This means 
that |?7m(l^)| < Cm with Cm a positive constant, on and inside the path C 
of integration in (j2.5p . Hence, \Um{\/t/n)\ < Cm on C (indeed, the path C 
has been modified after the change of variable, but we can set it equal to 
the previous path). Also, (1 — t/n)^^'^^^^ is bounded on C (and converges 
to e* for large n). Therefore, the above integral is 0(1) as n — t- oo. Hence, 

^™(^) = ^ (;^fe+T)7^ j ' ^ ^ oo, z / 0, (3.9) 

(n) 

which is comparable with the large n behavior of ^m ■ 

Prom the integral in (j3.6p it easily follows that <l?2A-+i ~ ^'-'^ ^ ~ 
0, 1, 2, . . . , n, and we see that, when n is large, only the even terms in the 

(n) 

series give contributions. We notice that constitute an asymptotic se- 
quence. This follows from 

We can collect the results of the section as follows. 
Theorem 1 For n ^ oo we have the asymptotic expansion 

oo 

Y>:{z) ~ {2zr 2^6^/^ {\-\k) , (3.11) 

fc=0 " 
which holds uniformly for \z\ > zq, where zq is a positive constant. 

In Table [1] we give the relative errors 5 when we use the expansion in 
()3.1ip with terms up and including k = 20, for // = 17/4 and several values 
of z and n. We see a quite uniform error with respect to z, except when 
z = ±^ 

Z ZC-^Q. 

Remark 1 If we wish we can expand the Pochhammer symbols in (13. IIP 
for large n and rearrange the series. In that way we can obtain an expansion 
of Yn{z) in negative powers of n, and this expansion is comparable with an 
expansion of which the first terms are given in (j3.ip . 



Remark 2 In (j3.1ip we expand the generalized Bessel polynomials Yn{z) 
in terms of another set of polynomials, the generalized Laguerre polyno- 
mials. Because the degree of these polynomials does not depend on the 
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Table 1: Relative errors 6 in the asymptotic expansion in (j3.1ip with terms 
up and including k = 20 for fj, = 17/4 and several values of z and n. 







z = 10^ 






z = -W 






j 


Y>:{z) 


6 


j 


Y^{z) 


5 


n = 50 


-1 


0.4232e034 


0.16e-03 


-1 


0.1961e026 


0.26e-ll 







0.1211e081 


0.38e-07 





0.1778e080 


0.62e-08 




1 


0.5131el30 


0.17e-07 


1 


0.4235el30 


0.14e-07 




2 


0.4707el80 


0.16e-07 


2 


0.4617el80 


0.15e-07 




3 


0.4666e230 


0.15e-07 


3 


0.4657e230 


0.15e-07 


n = 100 


-1 


0.1681e093 


0.30e-07 


-1 


0.5251e084 


0.68e-15 







0.3190el89 


O.lOe-10 





0.4501el88 


0.18e-ll 




1 


0.1325e289 


0.47e-ll 


1 


0.1089e289 


0.39e-ll 




2 


0.1213e389 


0.43e-ll 


2 


0.1189e389 


0.42e-ll 




3 


0.1202e489 


0.43e-ll 


3 


0.1200e489 


0.42e-ll 



large parameter, they can be evaluated much easier than the polynomials 
Yn{z). In fact, to compute the Laguerre polynomials we can use a recur- 
rence relation, which follows from differentiating ()3.3p with respect to W . 
Let Cfc = z), then 

z{k + l)ck+i = -{p.z + kz + l)ck - Ck-i, /c = 1,2,3, ... , (3.12) 

with initial values cq = 1, ci = —{^iz + 1)/ z. 

Remark 3 For general values of ^ the expansion in (j3.1ip is not conver- 
gent, but for ^ = 0,-1,-2,... it is. For example, for /i = a relation 
for the IT- Bessel function should arise. We have from 1^ L~''{l/z) = 
(— 1)'^/(A;! z*^), which gives the convergent expansion 

Y^iiz) = (2z)«eV^ E (l - • (3-13) 

A:=0 " 

Summing the series, separating the terms with even and odd k, we obtain 

Y^{z) = (-l)"eV^y^ (1/.) - /„^.(l/z)) , (3.14) 

and by using a well-known relation between the modified Bessel functions 
the representation in (12. ip arises. 
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Remark 4 The expansion in ()3.1ip is simpler than those of the following 
sections: it is easier to obtain and the coefficients are easily computed. When 
z is small the expansion in (j3.1ip breaks down. The Laguerre polynomials 
[k > 1) are not bounded, although the factor (2z)"' in front of the expansion 
has some control. But the main concern is the exponential factor e^/^, which 
has an essential singularity at 2: = 0. Recall that the polynomials Yn{z) all 
tend to unity when z — >■ 0. As mentioned in Theorem [H for the expansion 
in (j3.1ip we have to exclude a fixed neighborhood of the point z = 0. 



4 Expansions in terms of elementary functions 

By using saddle point methods we obtain expansions that hold uniformly 
inside sectors |ph it z| < ^vr — (5, where (5 is a fixed small positive number. 

For Yn{z) we take = as the large parameter. This gives a suitable 
identification of the results with those for the Bessel function K^.i (z) when 

/i = 0; see also (|2.ip . In addition we replace the argument z of Yn{z) by 
(observe that in [2] is also the large parameter, and the Bessel 
polynomial is considered with reversed argument). 

Because for both cases |ph ibz| < |vr — 5 the derivation of the asymptotic 
expansion is very similar we first summarize the results in the following two 
subsections, and in ^4.41 we give the details of the analysis. 

4.1 Expansion holding for |ph2| < ^tt 
Theorem 2 For large values of n we have the expansion 



(i+z2)i h 

and the expansion holds uniformly inside the sector jphzj < ^ir — 6. Here, 
6 is a small positive constant, v = n + ^, ( = l/{vz), Aq{^, z) = 1, 

A,i,..) = ^Z^-^^^^M^. (4.2) 



and 



_, (4.3) 
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For /i = the coefficients Ak{fi,z) reduce to those in the expansion in 
dZH), that is, Ak{0,z) = {-l)^Uklt), and 



>':(0 = f-^^-'K,A.z) ~ -^.---5: tlWQ, (44) 

which indeed gives the expansion in (j7.ip 

For C = ah Bessel polynomials Yn{() reduce to unity. We have as 

2: — ^ 00, t ^ 0, zt ^ 1, z - 7] ^ 0,. (4.5) 

As a consequence, Ai{fj,,z) — )• as C — ^ 0. In fact all coefficients Ak(fi,z) 
with k > 1 vanish as (" — )• 0, and both sides of ()4.ip reduce to unity. 

Recall that the simple expansion in ^is no longer valid when the argu- 
ment of the Bessel polynomials approaches the origin. 

4.2 Expansions holding for |ph(— z)| < |7r 

In this case we write (see (|2.12p and ()2.13p ) 

YJ^i-l/z) = FHil/z) + U^(\lz). (4.6) 
We have the following results. 
Theorem 3 For large values of n we have the expansions 



00 



1 + z + VTT^) 

k=0 



1 + z^) 4 ■ - 1^ 



(i + z- VTTz^Y \/z ^r, (ij 7) 

^ ' k=0 

and the expansions hold uniformly inside the sector |phz| < ^tt — 5. Here, 
Bo{fJ.,z) = 1, Co{fi,z) = 1, 

tjbt^ - 3) fit^z - 1) fi^jzt + l) 
B,{f,,z) = + , (4.9) 

and 

c,i,.o--'-^^'-^^^^- (4.10) 

The quantities t, and r] are as in Theorem\^ 

For ^ = the expansions reduce to those for the modified Bessel func- 
tions mentioned in (|2.14p . 
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4.3 Integral representations 

For deriving the asymptotic expansions we introduce the integrals 

/oo 
p^(s)e--<^W ds, 

(4.11) 



where > and 

(/)(s) =2zs-ln s-ln(l + s). (4.12) 

When 2: > the contour £ is a vertical line with 3fts > 0; when z is complex 
we can deform the contour in order to keep convergence. For the same 
purpose we can rotate the path of integration for Pv[z) in (14. lip . 

For certain choices of p^{s) and q^{s) these integrals give representations 
of the functions Yn[z), Fn{z), and Un{z). We have 

r(n + ^ + i) ^^^"^ - yi(rT^' ^^-'^^ 



(4.15) 



The multi- valued functions in cj){s), p^(s), and q^{s) have there principal 
branches and are real for s > 0. 

The representations in (I4.13P and (14.14p follow from the well-known in- 
tegral 

1 r°° 

U(a,c,z) = —- t''-\l + tY-''-^e-^^dt, Ka,z>0, (4.16) 

r(a) Jo 

the first line in (j2.1ip . and the second line in (|2.13p . For (j4.15p we refer to 
the first line in (|2.13p and the integral representation (see [101 P- 46]) 

1 ^ , . z^-^-e^ 



iFi(a; c; z) = ^— ^ / + s)"-'=s-» ds, (4.17) 



Tic) 

where £ is a vertical line with > 0. 
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4.4 Construction of the expansions 

We use the saddle point method for obtaining asymptotic expansions of the 
integrals in ()4.1ip . The saddle points follow from the equation (p'is) = 0, 
where 

^ (4,8) 

s[l + s) 

and are given by 



1- z + VTTz^ i-z-ViT 



2 



, = ■ . (4.19) 

When z > the saddle points are well-separated, with — 1 < s_ < — ^ and 
s+ > 0. We have the following limits: lim^^o s+ = +oo and lim2_j.oo s+ = 0. 
Also, 

s+(l + s+) = , (4.20) 

and 



l-z + lni2z) + v, ^"(s+) = -^^y==, (4.21) 

1 + V 1 + ^^'^ 



with T] defined in (|4.3p . 

We use Laplace's method with the transformation 

0(s) - 0(s+) = l^"{s+)w'^, sign{w) = sign(s - s+). (4.22) 
We have s = w + Oiw"^) as w — )• 0. The integrals in (|4.1ip become 



J-oo dw 



(4.23) 



Q'^(') = ^— 9He dw, 9iw)=q,{s) — 



By expanding /(w) = YlT=o fkW^ and g{w) = YlT=o dk'^'' obtain the 
asymptotic expansions 



Qi;(^)~5o^^W^;^^(-l) — ^ 



(4.24) 



fc=0 
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where (see also (|4.20p ) 



because ds/dw = 1 at it; = 0. 

By using (|il9]) - (|i:2T]) and (fiJi]) - (jiJ5]) it follows that 



(4.25) 



yA'fC)-^ — ^— - '■ \^Ek_Alhll (4 26) 



^^^^ ^ ' , .on! r*r„^,,^i^ 



j^k 



(4.27) 



/ . (i + z - vi + zA^ ^/^ ~ (7, r// z 



(4.28) 

The coefficients F^\^, z) are obtained from (j4.25p and the function f{w) of 

/OA 

(j4.23p with the function p^{s) as given in Km . and Fl '{ii,z) follow from 
taking the function Pfj,{s) as given in (j4.14p . The function T* is the slowly 
varying part of the corresponding gamma function. That is, 

r*i. + a)= f;^, 70(a) = 1, .^00. (4.29) 



k=0 



The coefficients 7/c(a) follow from standard methods for the gamma function; 
see S31 

The final form of the expansion of 1^(C) given in (j4.ip in Theorem [5] 
can be obtained by dividing the expansion in (j4.26p by the expansion of 
r*(i/ + // + ^) given in (I4.29p . and similar for the other expansions. This 
gives for A; = 0, 1, 2, . . . 

fc-i 

Ak{^l, z) = Fl:'\fi, z)-Y, Mf^^ (/^ + I) , (4.30) 

i=o 
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fe-1 

B.ifi, z) = z) - ^ z)7fc_, + I) , (4.31) 

i=o 

fc 

Cfc(^,z) = j;(-l)^G,(/i,z)7fe-, (i) . (4.32) 
i=o 

4.5 Computation of the coefficients 

To compute the coefficients -Ffe(/i, z) and Ak{n, z) we need the coefficients in 
the expansion 

oo 

s = s+ + ^SkW^, (4.33) 

k=l 

which follow from (j4.22p . We write, as in (j7.3p . 

t = , ^ (4.34) 



and obtain 



2-t (l-t)(5t3-6t2 + 2) 

si = 1, S2 = S3 - 



6 ' 18*2 
z{l - t)(40t^ - 65*3 _^ 24*2 _ 2t + 4) 
135*2 ' 



(4.35) 



With these coefficients we can compute the coefficients f{w) and g{w) of 
()4.23p by choosing the appropriate p^[s) and q^{s). 

To obtain the coefficients in (j4.30p - (|4.32p we first compute 7fc(/i + ^) 
that appear in (|4.29p . We have 

7o(m+|) = 1, 7i(m + |) = ^(-1 + 12^'), 

^2 = Tlk (1 + - 24^' - 19V + 144/) , 

73 = (1003 - 720^ - 17100/ + 11520/ + 32400/- 

34560/ + 8640/) , 

74 = 398]^ (-4027 - 288864^ + 151824/ + 1618560/- 

1239840/ - 1645056^^ + 2177280/ - 829440/i'^ + 103680/) . 

(4.36) 
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Figure 1: Saddle point contours of tlie first integral in (|4.1ip for z — e'^, 6 = fcvr/lO, 
k = 1,2,3,4,5. 




e= 5JI/10 



4.6 Extending the result to complex values of z 

From ^ p. 378] it follows that the expansions in ()7.ip and (j4.4p hold for large 
values of i' and are uniformly valid for complex values of z inside the sector 
Iphzl < ^vr — 5 with 5 a small positive number. As can be seen from the 
front factor and the coefficients, it becomes invalid when z approaches iti. 
In that case the singularities of the functions / and g in ()4.23p approach the 
origin. 

The singularities come from those of the mapping in (I4.22p . This map- 
ping does not depend on fj, and, hence, Laplace's method remains applicable 
for all fixed values of ^, and also (|4.ip is uniformly valid for complex values 
of z inside the sector |phz| < ^tt — 6. 

For complex z inside the sector jphzj < ^vr the saddle points given in 
(|4.19p move into the complex plane, and it is for all these values of z possible 
to find a single saddle point contour from through s+ such that ph(zs) = 
at infinity. If |phz| < ^vr — 5 the singular points of the mapping in (I4.22p 
and of the function / and g in (j4.23p stay away from the origin. 

In Figure [U we show the saddle points s+ (black balls) and the corre- 
sponding saddle point contours of the first integral in ()4.1ip for z = e*^, 
9 = /cvr/lO, k = 1,2, 3, 4, 5. The black balls at the left are the saddle points 
S-. When z = i the saddle points s+ and s_ coincide at s = + 
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The saddle point contours of the second integral in (j4.1ip are the paths of 
steepest ascent of the first integral. 



4.7 Comparison with earlier expansions 

In [2] the expansions of Yn{z) are given for large n with possibly large values 
of ji as well. This makes comparison with our expansions rather complicated. 
In fact in Dunster's expansions given in [21 §§6-7] the expansions can be 
re-expanded for small values of a parameter a corresponding with (in our 
notation) ^i/v, and in this way our results of the present section may be 
obtained. Dunster has used Olver's theory [8j for linear differential equations 
of second order, with bounds for the remainders in the expansions and a 
recursion formula for the coefficients. 

Starting from an integral we show how to include as a second large 
parameter, leaving out the details. Let = av and write (j5.4p in the form 



yM(^) = Sf!^ ^ / ^ ds, (4.37) 

where C = l/{i'z), u = n+ ^, and 

V'(s) = 2zs -{l + a)\ns- ln(l + s). (4.38) 

Then the saddle point analysis for the result of ^4.11 can be repeated, giving 
an expansion that holds again in the sector |phz| < ^ir — 6 and a > —1 + e, 
with 6, £ small positive numbers. In a similar way the expansions of §4.21 can 
be modified. 

With respect to the results in terms of elementary functions given in [13] 
we observe the following points. 

• A remarkable point is that for Yn{z) the point z = excluded, whereas 
in our results this point is accepted as long we approach it inside the 
sectors |ph(ib2;)| < ^ir — 6. In fact the results are essentially the same 
as our results, although the notation and scaling of the parameters is 
different. For example, the factor e~^/^ at the left-hand side of 
(2.26)] can be combined with e^""'"^^-^(''+'"^ to give a regular expression 
as z — 7- 0. Also, the domains of validity are different and are simpler 
in our case (just the sectors |ph(ib2:)| < hir — 6). 



• The expansions are derived from the integral representation (j2.5p , after 
a transformation. A detailed discussion is given about the location of 
saddle points and paths of steepest descent and the domains for the 
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expansions to be derived. Our start is from two integrals related to 
those for the Kummer functions. 

• The reader has to become familiar with several aspects of the detailed 
description of the paths of integration, the domains of validity and the 
proper choices of the branches of some multi- valued quantities. In this 
sense, our approach is more accessible for the reader who wants to use 
the results and to construct more terms. 

• Complex quantities arise in the coefficients and front factors (which, 
of course for real z will provide real expansions) . Our expansions show 
quantities that are real for real z. 

• It is not indicated how the results reduce to the well-known expansions 
of the modified Bessel functions; relations with the Kummer functions 
are not given. 

5 Expansions in terms of modified Bessel func- 
tions 

The expansions for l^(l/(z^2;)) in the previous section 21 become invalid 
when z approaches the points iti, because in that case the saddle points co- 
incide. As shown in [2j and |13j it is possible to derive uniform expansions in 
terms of Airy functions, and these expansions are valid in large z— domains. 

For the modified Bessel functions ly{yz^ and Kyiyz') similar asymptotic 
phenomena arise when z approaches the points ibi, and the expansion in 
terms of Airy functions is available in the literature. In fact expansions for 
the Hankel functions and the ordinary Bessel functions can be used. See |9l 
] and P Chapter 11]. 

Because the asymptotic phenomena of the generalized Bessel polynomi- 
als Yn{z^ for large n and fixed /i are the same as those of the polynomial 
we approach the problem for obtaining uniform expansions by ex- 
panding the generalized polynomials in terms of the modified Bessel func- 
tions Ky{^z^ (with V = n + ^), which are the same as the reduced Bessel 
polynomials Y^{z) (see (j2.ip ). By using the existing results for the Bessel 
functions a complete description is available in this way. 

We summarize the results of this section as follows. 
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Theorem 4 For n ^ oo we have the asymptotic expansion 



r(n + M+l)V » V ^ji-* ^jf'; 

and i/ie expansion holds uniformly with respect to all z. Here, C, = l/{vz), 
u = n+\, 

Co=f2-ie-^^V, Ci = :^(l-i)/i(M-l)(-2// + l + 3z)Co, (5.2) 
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and 

DQ = \{l-l)^xCo, D^ = -j-^i^i\^l-l){-^Ji + 2 + 2,i)Co. (5.3) 

5.1 The construction of the expansion 

To start the construction of the expansion we write (12. 4p in the form 
where again 



yM(^) = ^Af^^f^ / e-'^(-) d., (5.4) 



C = — , v = n + \, (/)(s) = 2zs - In s - ln(l + s). (5.5) 
vz ^ 

We write 

/o(s) = = ylo + B^s + 0'(5)5o(5), (5.6) 
and substitute s = and s = s_ to obtain 

s+/o(s-) - g-/o(g+) „ /o(g+) - /o(g-) 7X 

S+ — S- — S- 



Putting (j5.6p into (|5.4p we obtain 

KiO ^ A„.„ + B.., + «£M£) e-"*<.. (5.8) 

r(n + ^ + i)yo Vs(i + s) 

where 

.„ = ^?^v:((), ^. = '" v;'(c). (5.9) 

r(n + /i+l) r(n + ^+l) 
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By using ([2J]) and (pJO|) it follows that 



In the second line we can also write |1H p. 234] 

(1/z - l)K,{uz) - Kliuz) = K,+i{uz) - K,{uz), (5.11) 

but we prefer the notation with the derivative, because the asymptotic ex- 
pansions of Ku{vz) and Kl{vz) are quite related and usually presented to- 
gether. 

The next step is to use integration by parts in (|5.8|) . and this gives 

vT{n + + Jq ^s{l + s) 

where 

his) = V^(^)^^^7==v (^-13) 
ds ^s{l + s) 

Repeating this procedure by writing for A; > 

his) = Ak + BkS + Ms) = sf', (5.14) 

s+fkjs-) - s-fkis+) fk{s+) - fkjs-) , . 

Ak = , = , (5.15) 

s+ — S- •S+ — S- 

A« (») = V7(TT7)±-0= = ,t W - (5.10) 

we obtain for K > 

K-l . K~l 



k=0 fc=0 

(2z.z)"+/^+i r Ms) 



u^Fin + /i + 1) io ^^(^T7 



(5.17) 



We rearrange the expansion by using (j5.10p and writing 

1 - z 

to obtain the expansion given in Theorem [H 



Ck = Ak + ^^Bk, Dk = -\Bk, fc = 0,1,2,... (5.18) 
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Remark 5 To compute the coefficients , defined in (|5.15p and the 
functions say, by using a computer algebra package, it is convenient 

to write the functions fk{s) in the form of two-point Taylor expansions at 
the saddle points and S-. More details on this method can be found in 

Remark 6 For integer values of // we have the following simple cases. 

1. For /i = 0, 1, 2, . . . the expansion in ()5.17p has a finite number of terms 
which can also be obtained from the recursion in (j2.9p . 

2. For fi = —1,-2,— 3,... we can also obtain an exact result. When 
/i = — 1 we have 

Ck = ^, Dk = -^, A: = 0,1,2,..., (5.19) 

and we can sum the convergent series when 2z^ = 2n + 1 > 1. This 
gives a result that corresponds to the relation in ()2.9p with = —1. 

6 Concluding remarks 

In ^we have given a new simple expansion of Yn{z) that is valid outside 
a compact neighborhood of the origin in the z— plane and new forms of 
expansions in terms of elementary functions valid in the sectors |ph(ibz)| < 
^TT — (5 not containing the turning points z = -izi/n. To avoid mappings for 
obtaining expansions in terms of Airy functions we have given expansions in 
terms of modified Bessel functions. For these functions very detailed Airy- 
type expansions are available, which can be used to obtain similar expansions 

7 Appendix: Expansions of the modified Bessel 
functions 

Because we compare the expansions for the Bessel polynomials to those of 
the modified Bessel functions, we summarize a few details about the uniform 
expansions of the K— and /—Bessel functions. 
We have [sjl, p. 378] 

2http://dlmf.nist.gov/10.41 
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which hold when z/ — )• oo, uniformly with respect to z such that |phz| < 
^TT — 6, 6 being an arbitrary positive number in (0, ^tt). Here, 

The first coefficients Uk{t) are 

3t-5t^ 81*2-462*4 + 385*6 
no(*) = l, ^i(t) = , n2(t) = ^^^^ , (7.4) 

and other coefficients can be obtained by applying the formula 

Uk+i{t) = 1*2(1 - *2)4(t) + I /" (1 - 5s^)uk{s)ds, A; = 0, 1, 2, ... . (7.5) 

Jo 

For the derivatives we have 

* fe=o 

where 

, -9* + 7*3 -135*2 + 594*4-455**5 
vo(i) = 1, vi{t) = — , V2{t) = , 

(7.8) 

and other coefficients can be obtained by applying the formula 

t'fc(t) = nfc(*)+*(*2-l)(i7Xfe_i(*) + *4_i(*)), A: = 0,1,2,.... (7.9) 

The expansions in (|7.ip and (|7.2p become invalid when z approaches the 
turning points ±i. In that case expansions are available in terms of Airy 
functions. First the functions luiz) and K^{z) should be written in terms 
of ordinary Bessel functions, and then the results for these functions can be 
used; see [H p. 419-426]. 
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